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Abstract 



"^ . The summed strength for transitions from the ground state of ^^C 

[JL( ! via the operators st, it, rY't, r\Y' s\^t and r[y'£]''^t are calculated using 

the AA^ = rotational model. If we choose the z component of the 
isospin operator t^, the above operators are relevant to electromag- 
netic transitions; if we choose i+ they are relevant to weak transitions 
^^ I such as neutrino capture. In going from the spherical limit to the 

asymptotic (oblate) limit the strength for the operator st decreases 

steadily to zero; the strength for the operator £r (scissors mode) in- 

(^ \ creases by a factor of three. For the last three operators - isovector 

^ ' dipole, spin dipole and orbital dipole (including the twist mode) it 

is shown that the summed strength is independant of deformation. 
The main difference in the behavior is that for the first two operators 
we have in-shell transitions whereas for the last three operators the 
y . transitions are out of shell. 
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1) Introduction 

The present generation of neutrino-oscillation experiments [1] require the 
knowledge of neturino-nucleus reaction cross-sections of both exclusive and 
inclusive processes. In particular, neutrino reactions on the ^^C nucleus have 
been very instrumental in studying neutrino oscillations in the LSND [2] and 
KARMEN [3] experiments. 

Neutrino flux averaged ^'^C{i'^, fi^) and ^'^C{i'e, e") inclusive cross section 
as well as the exclusive ^"^Ci^u^, ji'Y'^Ngs. and ^^C(z/e, e~)^^A^c,.s. were deter- 
mined in the above experiments. Theoretical and experimental attention has 
been given to the ^^C(z/^, /i~) measurement performed at LSND with a flux of 
neutrinos obtained from the decay in flight (DIF) of pions. The experimental 
flux averaged cross-section to the ^^A^ ground state is [4]. 

ae^e= (0.66 ±0.1 ±0.1) X lO'^^cm^ (1) 

while the flux averaged total inclusive cross-section (which includes the tran- 
sition to ^^A^g.s. is now put at [4], 

ainc = (12.2 ± 0.3 ± 1.8) x lO^^^cm^ (2) 

When the ^^C ground state is assumed to form apg^ closed shell the Hartree- 
Fock (HF) and the subsequent Random Phase Approximation (RPA) calcu- 
lation give an exclusive cross-section which is a factor 5-7 [5,6] too large 
compared to the LSND result. The inclusive cross-section in the RPA turns 
out to be (20 — 22) x 10~^^cm^ [5,6] almost a factor of two larger than the 
LSND result. However, it is well known for a long time that the ^^C ground 
state is not well described by a pure pl,2 conflguration. The Gamow- Teller 
transition to the J = 1'^ g.s. of ""^^A^ as well as the Ml transition in ^^C 
are too large by factors of about 5 when a pure p^^ conflguration is used 
for the ^^C g.s. Conflguration mixing and the use of intermediate coupling 
wave functions reduces these transitions. In the early theoretical work [5] the 
authors have accounted for the fact that the Pg ^ g.s. is not a good starting 
point by simply dividing the calculated exclusive {u^ji^) cross section by 
approximately a factor of flve. They have however not altered the inclusive 
cross-section to the excited states, leaving unchanged the RPA result, which 
was calculated starting from a pure pg /2 g.s. In ref. [6] the reduction of the ex- 
clusive cross-section was achieved by using a pairing model approximation. 



To account for the configuration mixing in the g.s. of ^^C a Bogohubov- 
Valatin approximation was apphed and partial occuppation numbers for the 
P3/2 and pi/2 orbits were assigned, {vp.^,^ = 1, Vp^.^ = 0). This leads to a large 
reduction in the exclusive cross section which is determined completely by 
the Gamow- Teller strength of the J = 1'^ g.s. in ^^A^. In a very approximate 
way the effects of pairing on the inclusive cross-sections were also estimated. 
These estimates, which included also the energy shift due to pairing, led to 
a 20-25% reduction. 

In this work we wish to study the effects of deformation on the transitions 
from the ground state of ^^C Indeed ^^C is generally regarded as being 
strongly deformed oblate nucleus and a deformation parameter (3 ~ —0.6 
is often quoted. One cannot simulate deformation effects by doing an RPA 
calculation on a spherical ground state. In nuclei where pairing correlations 
are important, both deformation and pairing should be taken into account. 
It is not clear how important pairing is for the light nuclei. 

We feel that not too much is known on the effects of deformation on 
isovector transitions - the main focus has been on the isoscalar enhancement 
of the B(E2) Oi -^ 2, in ground state rotational bands. We therefore shall 
focus our attention on simple isovector operators which are relevant to both 
electromagnetic and weak transitions, and we shall consider the simplest 
problem of summed strengths to states (A, K) where A is the final total 
angular momentum and K the projection of total angular momentum on the 
symmetry axis in the intrinsic state of a deformed rotating nucleus. 

2) ^^C in the Spherical and Asymptotic Limits 

In the spherical limit ^^C consists of 4 nucleons in the Osi shell and eight 
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in the ps shell. An elastic dipole transition could take place either by exciting 
a nucleus from opa to the Is-Od shell or by exciting a Osi nucleon to P3/2. 
An Ml transition would consist of exciting a ps/2 nucleon to the pi/2 shell. 

It was shown very early on e.g. by Cohen and Kurath [7] that this limit 
is a very bad approximation. There is considerable configuration mixing. In 
the spherical limit the J = 1"^ state wuld consist of seven P3/2 and one pi/2 
nucleon. However, in shell model calculations one finds that the number of 
Pi/2 nucleons in the J = if state is almost the same as in the ground state. 

We next consider the asymptotic oblate limit. It is generally believed 
that the ground state band in ^^C is oblate with f] ~ —0.6. Indeed in the 
deformed oscillator model with the Mottelson conditions one obtains such a 



result. We may therefore expect that this hmit is closer to reality than the 
spherical limit. 

In the asymptotic oblate limit we have the following single-particle levels 
and wavefunctions. 
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{K = \) 


occupied 


2) 
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rFiiT 


m-i i) 


(^ = 1) 


occupied 


3) 


N=l 
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m-i T) 


K = \ 


occupied 


4) 


N=l 


rVioT 


(^0 i) 


K = \ 


empty 



We can see immediately that the matrix element of the spin operators s 
or st from the J = 0^ ground state to an excited J = 1+ state in ^^C will 
vanish. The spin operator cannot change Yn into Yiq. The result is more 
general-in the SU(3) limit the spin Mi's will vanish. 

3) The Ml Transition in ^^C 

We quote the results and use the conventions of reduced motion elements 
of Bohr and Mottelson [8] . 

Ml = i/^/2 (3) 



<J2 = £+^||Ml||ji = £-l> 




^^(^.-(?,)^2ji + l(ji-10|j22) (4) 

Thus < p± 1 1 Ml I |pi >= — ^ (s's — Qf)- Here the bare values of the parameters 
are: gi^ = l,geu = 0,gsn = 5.586, gsu = -3.826. 

For particle-hole matrix elements of a tensor operator the same reference 
[8] gives: 

In the spherical limit the J = 1+,T = 1 state in ^^C is 
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Our result is then 

B{Ml) = j-—[{gsn - 9sv) - {gin - giv)f 

(27r) 
= 11.262/i^ (7) 

In the rotational model of Bohr and Mottelson the following result is 
given for an Ml transition from a J=0, K=0 state 

B{Ml) = \{J = 1+K\M1\J = 0K = 0)p(2 - 6k,o) (8) 

the factor of 2 for K not equal to zero comes from the fact that one can excite 
both K and -K. 

For an intermediate deformation the levels of the previous section become: 

2) N = 1 

3) N = 1 

4) N = 1 



Fnt 


occupied 


cFio T + dYu i 


occupied 


-dY,o T + cYn i 


empty 



In the asymptotic oblate limit c = 0,d = 1. In the spherical limit c 
One gets both K=l and K=0 contributions to Ml 



K = 1 a) <4\Ml\3>=-cdg'^ + {c^ -d'^)g'LV2 

K = I h) <4:\Ml\2>=cg'^-V2dg'^ 

K = c) < A\Ml\3 >= -cd{g'^ - g'^) 

K = d) <l\Ml\3>=cd{g'^-g'L) 



Here 



(fi'sTT — gsu) I _ ge-TT — giu 
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d). 



V2 "^ v^ 

To get the total B(M1) we add up the sum of the squares of a) b) c) and 



B{Ml) = ^n{[-cdg', + ic'-d')g',V2r 

+ [cg:-V2dg',r + 2c'd'ig:-g',f} (9) 

In the above expression the first two terms come from excitations to K=l 
states and the last to K=0 states. This becomes identical to the spherical 
expression previously given for c = w |, dJ^. 
In the asymptotic oblate limit this becomes 

B{M1) = ^{gL^ - gL.f = 0.4775/i2, (10) 

ZTT 

Only the orbital part of Ml contributes in this limit. This is known as a 
scissors mode excitation. 

We also consider an intermediate case with 5 = —0.3 where we use the 
wavefunctions of Max Irvine [9]. For these c=0. 25457 and d = 0.96706 we 
obtain B(M1) = 2.655 fi% which is close to the experimental value B{M1) = 



2.61/x? 
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We now consider the operator relevant to the neutrino capture z/^+ C — >■ 
N + fi~ . For the Ml allowed we consider the operator 2st+. We can get 
the result for this from the previous section by dropping the facrtor jTt and 



by setting g'^ to 2 and g'f^ to zero. Hence B{S = 1, L = 0) = I2c^d'^ + 4c^. In 
the asymptotic oblate limit c=0 so this vanishes. In the spherical limit this 
becomes y. For an intermediate deformation, 6 = —0.3 the value is 0.987 
(we do not include a coupling constant Cgt = 1.251 in our definition of B). 

4) The operators rYi^fct+,r[Fis|^t+ and r[Yi£]^t+ 

At finite momentum transfer we can get other contributions to the neu- 
trino cross section. We here consider the operators listed in the title. The 
operator rYi^^t^ arises from the "Forbidden Fermi", r[Yis]^t+ from "forbid- 
den Gamow- Teller" . The dipole orbital operator is included for completeness. 

We give results for the total strength for these operators in Table 1. We 
will give less details than in the last section. The following formulae are 
useful for constructing the table. 

a) Expansion of the Dipole Spin Operator. 

[Y\s]: = -l/VSY^So + 1/VQYIs^ - 1/VqY1,s+ 



[Y\s]l = ir^_ + ir\s+ 

V2[Y's]\ = Ylso + ^Y}s^ 
V2[Yh]l = Yls^ - ^Y^s^ 
V2[Yh]l = -Y,Ks+ 

Here 

<TSoT>=| <isoi>=-| 

<t S+ i>= 1 <i s_ t>= 1 

b) The following radial integral is useful 



/■ y^i y^2 y^3 ^, . _ r (2^i + 1)(2^2 + 1)(24 + 1) ni ^e^e^e-,) liii^h 

/ mim2 ms^^ [ a J \o o o / Vm,im2'' 
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where the last two factors are 3j symbols. The following radial integrals are 
needed 

J OS r op r'^dr = J^b 

J op r od r'^dr = J^b 
J op r Is r^dr = b 

where harmonic oscillator wave functions have been used and b = Jh/mto. 

As in the previous section we consider transition to states with quantum 
numbers J and K. We form all particle-hole states 

Kp,KHiK = Kp-KH)K>0 

We then calculate the matrix element M{Kp, Kh, J)k- For given K, J 
we get the total strength S(J,K) such that 

S{J,K)= Y. \M{Kp,Kh,J)\' (12) 

KpKh 
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Note that we consider not only excitations from Op to Is ,0d but also from 
Os to Op. To show the results in as simple a way as possible we multiply the 
strengths (i.e. squared matrix elements) by Anmuj/h. 

The results, which are shown in Table 1 are very simple and very striking. 
For the operators rY^t, r[Y^s]^ and r[y^£]'^ the total summed strengths 
are exactly the same in the spherical limit and the asymptotic limit. This 
is quite different from the operator st, the strength of which is large in the 
spherical limit but vanishes in the asymptotic limit. Also for the operator ir 
(scissors mode) the strength in the asymptotic limit is different then in the 
spherical case - it is a factor of three larger. 

Note that for a given A in the spin-dipole case [rF^s]^ the summed 
strength is not the same in the two limits. In the asymptotic case the partial 
sums for A = 0, 1 and 2 are proportional to (2A + 1) - 2.25, 6.75 and 11.25. 
In the spherical case they are 3.25, 8.25 and 8.75 respectively. But the total 
sum is the same in the two cases. 

The redistribution of the dipole strength (spin independent and depen- 
dent) and energy shifts arising from deformation may lead to some modifi- 
cations in the inclusive neutrino cross-sections. This is due to the fact that 
there is energy dependence in the neutrino fluxes used in the experiment and 
also due to the energy dependence of the neutrino-nucleus interaction. The 
basic conclusion that the dipole strength is not affected by deformation is 
a very strong indication that the deformation will not change the inclusive 
cross-section a lot while at the same time the exclusive cross-sections to the 
^^A^ g.s. will be reduced considerably by the deformation. 

The Cal-Tech group [5] did not include deformation effects explicitly but 
simulated the consequences of these decreasing the "allowed" GT strength 
in the neutrino raction by a factor of five without altering the total dipole 
strength. Our calculation justifies this in part but there are residual effects 
of deformation on the dipole strengths in the sense that the strengths to 
indivdual J values do depend on deformation. Since the various J components 
are split in energy this could affect the neutrino cross section. 

Note that in the dipole-orbital case there is no contribution for A = 0. The 
operator for this case can be written as Cf- i where C is a constant. We can 
write f-i = f-fxp = {fxf)-p. This vanishes because fxf= 0. Alternatively 
we can say that r ■ i given us the component of angular momentum in the 
radial direction - this is zero. 



Note that in the case [rY^i]^^'^t one excites the twist mode. Like the other 
cases we get the same results for this mode in the spherical and asymptotic 
limits. 

5) Proof of the Independence of the Summed Strengths for the Operators 
j^yi^ j^jyi^^A g^j^^ [r'K^f]''* on Deformation. 

We consider the most difficult case the spin-dipole. We consider transition 
from a A; = I orbit xY^ | +yYi |. The summer strength is 

S= Yl \<LMMs[Y's]Uyo^+yyll>\' (13) 

LMMsXK 

We expand [Y\s]^^ = E„;3(ll«/3|Air)y>^. 

We use the orthonormality property of the Clebsh-Gordan coefficients 

J2{ilaP\Xk){na'\P'Xk) = 6^^.6pp. 

Xk 

to obtain 

S= Y. \< LMM.Y^spxY^ t +yyl i> r 

We can evaluate the spin part using <\ So\>= \ <[ So [>= —\ <] 



S = S{Ms = ^,(3 = 0) + S{Ms = ^,f3 = l) 

S = Y{i\<LMY^yo>\' + ^\<LMY^yi>\' 

aLM ^ ^ 

+ jl < LMYlYl > I' + yl < LMY^Yl > \^} 



By using the Wigner-Eckart theorem we can show 
Y I < LMY^Yf^ > |2 = ^ I < LMY^Yi > p = i| < L\\Y'^\\Y^ > 

Ma Ma "^ 
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Hence the summed strength becomes 

^ L 

However for the orbit x Y} f +yYi | the normahzation condition is x'^+y'^ = 
1. Hense S is independent of x and y and hence of deformation. 

Note that the independence of deformation holds also separately for each 
L. The proofs for the other modes rV^ and r[y^£]'^ go along the same lines 
as above and are in fact much simpler. 

6) Summed Strength - An Alternate Approach 

We now show that the summed strength is independent of the shell model 
configuration as long as all valence particles in the ground state are in the 
Op shell and that the expectation value of r.r is the same for all particles in 
the p shell. 

Let the dipole operator be O. The summed strength is the ground state 
expectation value of Sum(i,j)0(i).0(j). If i does not equal j we get zero if all 
nucleons are in the Op shell because of parity. If i equals j then if 0=RY^ 
we get basically r.r and this is the same for all particles. For the spin dipole 
we get for J=0 s.r s.r which equals 3/4 r.r. For J=l we get s ■ f s ■ r which 
equals 1/2 r.r. For J=2 we get (ss)J=2.(rr)J=2. But consider (ss)J=2 Ms=2. 
This equals 1/2 s+s+. This vanishes because we cannot raise the spin of a 
spin 1/2 particle twice. The other components of (ss)J=2 will also vanish. 
Thus the ground state expectation value is reduced to one body terms of the 
form r.r which will be the same for P3/2 and pi/2 particles provided harmonic 
oscillator wave functions are used. To get model dependent results we have 
to include admixtures from higher shells and/or use non harmonic oscillator 
radial integrals. 

To make the argument complete we must also take into account the fact 
that there are particles in the Os shell and one can have matrix elements of 
0(i).0(j) between Os and Op. This will have the structure of an exchange 
single particle interaction of a Op particle with a closed Os shell. This single 
particle matrix element will be the same for ^3/2 as for pi/2 only a 2 body 
spin-orbit interaction can lead to a single particle splitting of ^3/2 and pi/2- 

Extensive experimental work on ^^C(n, pY'^B reaction has been performed 
by N. Olsson et. al [10]. Recently a preprint has appeared in the LANL 
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network by P. Vogel [11] in which detailed calculations show that the total 
strength does not depend on pi/2 shell occupancy. 

This work was supported by the U.S. Department of Energy under grant 
DE-FG02-95ER-40940. 
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Table 1 Isovector dipole strength in ^^C{^^B{X)0+ -^ A) 

Dipole rY^t 

Asymptotic Spherical 

9 9 

18 18 



A 


K 








1 


1 


Sum 



27 27 



Spin Dipole 


r[Y^s]H 




A K 









2.25 


3.25 


1 


2.25 


2.75 


1 1 


4.50 


5.50 


2 


2.25 


1.75 


2 1 


4.50 


3.50 


2 2 


4.50 


3.50 


Sum 


20.25 


20.25 



Orbital 


Dipole r[YH]H 




A K 















1 


3 


6 


1 1 


11 


8 


2 


9 


6 


2 1 


19 


16 


2 2 


6 


12 



Sum 48 48 
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Table 2 Same as Table 1 but only for Os -^ Op Transitions 

Asymptotic Spherical 



0.75 
0.50 
1.00 
0.00 
0.00 
0.00 



A K 


[rY'S]'^ 





0.25 


1 





1 1 


0.75 


2 


0.5 


2 1 


0.75 


2 2 


0.00 



rY^ 



X K 




1 


3 


1 1 






2.25 2.25 



1 

2 
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